Ab-initio study of the electric transport in gold nanocontacts 
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By employing a real-space formulation of the Kubo-Greenwood equation based on a Green's 
function embedding technique combined with the fully relativistic spin-polarized Korringa-Kohn- 
Rostoker method a detailed investigation of the electrical transport through atomic-scaled contacts 
between two Au(001) semi-infinite systems is presented. Following a careful numerical test of the 
method the conductance of Au nanocontacts with different geometries is calculated. In particular, 
for a contact formed by a linear chain of Au atoms a conductance near 1 Go is obtained. The 
influence of transition metal impurities (Pd, Fe and Co) placed on various positions near the center 
of a particular contact is also studied. We found that the conductance is very sensitive to the 
position of the magnetic impurities and that the mechanism for the occurring relative changes can 
mainly be attributed to the impurities' minority d-band inducing resonant line-shapes in the s-like 
DOS at the center of the contact. 
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I. INTRODUCTION 



i The number of theoretical and experimental investigations of electronic structure and transport properties of atomic- 
sized conductors has greatly been increased over the last decade. The increasing interest for investigating atomic- 
sized conductors is driven by the possibility of using such systems in future nanoelectronic technologies. Widely 
applied methods for fabricating nanocontacts between macroscopic electrodes are the mechanically controllable break 
junction (MCBJ) techniquei* 2 *^ and scanning tunneling microscopy (STM)£*£*£ by pushing the tip intentionally into 
the surface. The crucial problems for both methods are the presence of contaminants and the mechanical stability. 
At sufficiently low temperatures the measurements revealed a quantized conductance for atomic sized nanocontacts 
made of various materials, not only pure metals but also alloys^. Nanocontacts made of gold are presumably 
the most studied systems in the literature both theoretically and experimentally. A dominant peak very close to 
the conductance quantum, 1 Go — 2e 2 //i, has been reported for gold (and other noble metals) in the conductance 
histograms*!, attributed to the highly transmitting sp-channel across a linear chain connecting the two electrodes. 
It was also found that the chain formation is in close connection with surface reconstruction phenomenal. For a 
comprehensive review of the field of atomic-sized conductors, see Ref. 

In order to understand the mechanism of nanocontact formation, electronic structure and transport, different theo- 
retical methods have been developed. Some theoretical studies use tight-binding method a 11 ! 12 , others are based on 
ab initio density functional theorjii&iii 5 .. Most of the transport studies rely on the Landauer-Biittiker approacbiSiii, 
although Baranger and Stone adopted the more sophisticated Kubo-Greenwood formulai&i2i22*2i for calculating the 
conductance between free electron leads^S,. By using this approach a recent study— focused on the effect of transition 
metal imperfections inserted into an infinite Cu wire showing that the conductance of the wire decreased due to the 
different conductance for the two spin channels (spin- filter effect). 

The fully relativistic screened Korringa-Kohn-Rostoker (SKKR) Green's function method proved to be an effective 
method to calculate electronic structure and magnetic properties of layered systems^iSi. Combined with an em- 
bedding technique based on multiple scattering theory, calculations have been performed for magnetic clusters on 
surfaces 25 - 26 27 . Employing the Kubo-Greenwood formula within this method permits one to investigate transport 
properties of atomic sized structures 2 ^. In this paper we report on calculations of transport properties of gold nanocon- 
tacts. We first briefly review the theoretical background of the applied method and give numerical evidence of its 
reliability. We then calculate and analyze the conductance for different atomic arrangements between semi-infinite 
Au(001) systems and investigate the effect of transition metal impurities on the conductance. We find a qualitatively 
satisfactory explanation of the observed changes in the conductance in terms of changes in the s-like local density of 
states (LDOS) at the center of the point contact caused by interactions with the d-like states of the impurity. 
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II. THEORETICAL APPROACH 



A. Expression of the conductivity 



The static limit of the optical conductivity tensor is given by the Kubo-Luttinger formulai&Si as 

h 
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where V is the volume of the system, /(e) is the Fermi-Dirac distribution function, J M is the iiXh component (p = x,y 
or z) of the current density operator, G ± (e) are the corresponding (upper or lower) side limits of the resolvent of an 
appropriate Kohn-Sham(-Dirac) Hamiltonian, while Tr denotes the trace of an operator. Integration by parts yields 

de^-S^ie), (2) 



with 



5 ^ = -2^/_>' 

x Tr | jj^lj v [^) - G~( e ')] - UG + {e') CT(e')]J l/ ^^ , (3) 

which has the meaning of a zero-temperature, energy dependent conductivity. For T = 0, <j^ u is obviously given by 

0>„ = (e F ) ■ (4) 
A numerically tractable formula can be obtained only for the diagonal elements of the conductivity tensor, namely, 

Tr{j,^±J,[G+(s)-G-(e)] J»[G+(e) - G~(e)] ^^^} = 
= Tr | [G+(e) - G~( e )] J M [G+(e) - G~( e )]} 

= ~Tr{J /1 [GH-( £ )-G-( e )]J„[G+(e)-G-( E )]} , (5) 
yielding the widely used Kubo-Greenwood formul a 19 ' 21 of the dc-conductivity at finite temperatures, 



On the other hand, Eq. (JTJ can be reformulated as follows 
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namely in terms of an equation which is similar to the formulation of Baranger and Stoned but clearly can be cast 
into a relativistic form. 
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B. Expression of the conductance 



Linear response theory applies to an arbitrary choice for the perturbating electric field because the response function 
is obtained in the zero limit of perturbation. Let us assume, therefore, that a constant electric field, E^, pointing 
along the z axis, i.e., normal to the planes, is applied in all cells of layer J. Denoting the z component of current 
density averaged over cell i in layer / by j l J , the microscopic Ohm's law reads as 

£ = ^-'E°-% iJ El, (8) 

Vat 

J 

where V a t is the volume of the unit cell in layer I. Note, that in neglecting lattice relaxations, V a t is uniform in the 
whole system. According to the Kubo-Greenwood formula, Eq. the zz component of the non-local conductivity 
tensor, can be written at zero temperature as 

4i' jJ = -^[ d\ u f d\' 3j (9) 
x Tr (j z [G+(e F ; v u , r^) - G~(e F - v u , rj v )] J z [G+(e F ; i> v , v u ) - G~ (e F ; r'^vu)]) . 

Here the integration is carried out over the ith unit cell in layer J, flu, and the jth unit cell in layer J, Cljj, while 
Tr denotes a trace over four-component spinors. The total current flowing through layer / can be written as 

Itot = A^j" = gll , (10) 

i 

where the applied voltage U is 

U = E J z d x , (11) 

and An and d± denote the area of the 2D unit cell and the interlayer spacing, respectively (V a t — A\\ d±). Combining 
Eqs. (|8|). I|1U|) and (|ll|l results in an expression for the conductance, 

s = ^££^ iJ . as) 

- 1 - i j 

where the summations should, in principle, be carried out over all the cells in layers I and J. An alternative choice of 
the non-local conductivity tensor is given by Eq. J7J). This leads to a huge simplification for the conductance because, 
as shown by Baranger and Stoned for free electron leads, the second term appearing in Eq. J7J becomes identically 
zero when integrated over the layers, I ^ J. It should be noted that recently Mavropoulos et al. 29 rederived this 
result by assuming Bloch boundary conditions for the leads. The conductance can thus be written as 
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2-Kd 2 , 



EE / rf3 ^ / d^jTr^G+ie^vuy^G-ieF^j.vu)] . (13) 



It has to be emphasized that because of the use of linear response theory and current conservation, the choice of layers 
I and J is arbitrary in the above formula. The numerical test of the method will clearly demonstrate this feature 
(see Section III.). On the other hand, as shown in R.ef. I29l when the layers / and J are asymptotically far from each 
other, the present formalism naturally recovers the Landauer-Biittiker approachiS^i. 



C. Computational details 

Using the embedding technique of multiple scattering the matrix representation of the scattering path operator (SPO) 
of a given cluster, Z clus { e ) — = | t qq' ( e )}> w i tn * an d 3 denoting sites in the cluster, Q and Q' indexing 

angular momentum quantum numbers, can be expressed as2^ 

(14) 
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where t (e). t (e) stand for the corresponding single-site t-matrices of the host medium and the cluster, and 

=host K 11 =clus y J 1 to to ' 

r, (e) for the host-SPO. For a two-dimensional (2D) translational invariant host medium, the latter one is calculated 

=host v 7 x 7 

by 

(e) = ^-f d^e^-^riUe^) , (15) 
"sz Jbz 

where Tj and Tj are 2D lattice vectors and the integral is performed over the 2D Brillouin zone of area flsz- 

The self-consistent calculations for both the host and the finite clusters were performed within the local spin-density 
approximation (LSDA)— , by using the atomic sphere approximation (ASA) and l m ax = 2 for the angular momentum 
expansion. The semi-infinite host system was evaluated in terms of the screened Korringa-Kohn-Rostoker method 
(SKKR)2i24 by sampling 45 fcii points in the irreducible (1/8) part of the fcc(OOl) Brillouin zone, see Eq. Q15|). and 
16 energy points along a semi-circular contour in the upper complex energy semi-plane. The latter set-up also applied 
for the self-consistent cluster calculations, whereby a sufficiently large number of atoms from the neighboring host 
(including sites in the vacuum) was taken into account in order to serve as a buffer for charge fluctuations. In the case 
of magnetic impurities, the orientation of magnetization was chosen to be normal to the fcc(OOl) planes (direction z). 
Additional calculations of the magnetic anisotropy energy confirmed this choice. 

In terms of the SPO the conductance in Eq. I|13fl can be calculated as 

-L i j 

where J^(e,e') stands for the relativistic angular momentum representation of the current density operator in cell i 
of layer / (see, e.g., in Ref. 28), and, correspondingly, the trace is performed in angular momentum space. Inherent 
to the SKKR method, a finite imaginary part, 5, of the Fermi energy has to be applied, = £p ±iS, which, however, 
has to be continued to zero in order to ensure current conservation. Concomitantly, the number of k\\ points taken 
in Eq. I|15(l has to be considerably increased. All results presented in the next Section for the conductance refer to 
5 = 1/iRyd. 

In the present work no geometry optimization has been carried out, that means all of the considered sites (both Au, 
vacuum and impurity sites) correspond to the positions of an ideal fcc(OOl) structure of gold with a lattice constant 
of asp — 7.68 a.u. A schematic view of a typical contact is displayed in Fig. ^ As follows from the above, atomic 
sites refer to layers for which we use the notation: C the central layer, C — 1 and C + 1 the layers below and above, 
etc. For the contact shown in Fig. e.g., the central layer contains 1 Au atom (the rest is built up from empty 
spheres), layers C — 1 and C + 1 contain 4 Au atoms, layers C — 2 and C + 2 contain 9 Au atoms and, though not 
shown, all layers C — n and C + n (n > 3) are completely filled with Au atoms and will be denoted by full layers. 



III. RESULTS AND DISCUSSION 



A. Numerical tests on different gold contacts 



As mentioned in Section II a finite Fermi level broadening, 5, has to be used for the conductance calculations. As 
an example, for the point contact depicted in Fig. we investigated the dependence of the conductance on S. The 
summation in Eq. i|16|) was carried out up to convergence for the first two (symmetric) full layers (/ = C— 3, J = C+3). 
As can be seen from Fig. [31 the calculated conductances depend strongly but nearly linear on 5. A straight line fitted 
for S > 1.5 mRyd intersects the vertical axis at 2.38 Go- Assuringly enough, a calculation with 6 = l^Ryd resulted in 
g = 2.40Go- Although the nearly linear dependence of the conductance with respect to 5 enables an easy extrapolation 
to 5 = 0, as what follows all the calculations refer to S = 1/iRyd. 

For the same type of contact we investigated the convergence of the summation in Eq. I|16|l over the layers / and 
J, whereby we chose different symmetric pairs of full layers. The convergence with respect to the number of atoms 
in the layers is shown in Fig. Convergence for about 20 atoms can be obtained for the first two full layers 
(/ = C — 3, J = C + 3), whereas the number of sites needed to get convergent sums gradually increases if one takes 
layers farther away from the center of the contact. This kind of convergence property is qualitatively understandable 
since the current flows from the contact within a cone of some opening angle that cuts out sheets of increasing area 
from the corresponding layers. As all the calculations were performed with 5 = 1/iRyd, current conservation has to 
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be expected. Consequently the calculated conductance ought to be independent with respect to the layers chosen for 
the summation in Eq. I|ltj|) . As can be seen from Fig. 0] this is satisfied within a relative error of less then 10%. It 
should be noted, however, that for the pairs of layers, I — C — n, J — C + n, n > 6 convergence was not achieved 
within this accuracy: by taking more sites in the summations even a better coincidence of the calculated conductance 
values for different pairs of layers can be expected. Fig. 21 also implies that an application of the Landauer-Biittiker 
approach to calculate the conductance of nanocontacts is numerically more tedious than the present one, since, in 
principle, two layers situated infinitely far from each other have to be taken in order to represent the leads. 

Although only one Au atom is placed in the center of the point contact considered above, see Fig. the calculated 
conductance is more than twice as large as the conductance unit. This is easy to understand since the planes C — 1 
and C + 1, each containing four Au atoms, are relatively close to each other and, therefore, tunneling contributes 
quite a lot to the conductance through the contact. In order to obtain a conductance around 1 Go, detected in the 
experiments, a linear chain has to be considered. The existence of such linear chains is obvious from the long plateau 
of the corresponding conductance trace with respect to the piezo voltage in the break-junction experiments. Since, as 
mentioned in Section II, our method at present can only handle geometrical structures confined to three dimensional 
translational invariant simple bulk parent lattices, as the simplest model of such a contact we considered a slanted 
linear chain as shown in Fig. ^jp. In here, the middle layer (G) and the adjacent layers (G ± 1) contain only one Au 
atom, layers G ± 2 and G ± 3 4 and 9 Au atoms, respectively, while layers G ± 4 refer to the first two full layers. 
The sum in Eq. Ijltil) was carried out for two pairs of layers, namely for I = C — 4, J = C + 4 (full layers) and for 
I = C — 2, J = C + 2 (not full layers). The convergence with respect to the number of atoms in the chosen layers can 
be seen from Fig. [S] The respective converged values are 1.10 Go and 1.17 Go- In the case of I = C — 2, J = G + 2 
we observed that the contribution coming from the vacuum sites is nearly zero: considering only 4 Au atoms in the 
summation already gave a value for the conductance close to the converged one. The small difference between the two 
calculated values, 0.07 Go, can most likely be attributed to an error caused by the ASA. Nevertheless, as expected, 
the calculated conductance is very near to the ideal value of 1 Go- 

Another interesting structure is the 2x2 chain described in Ref. 14 Here we considered a finite length of this 
structure sandwiched between two semi-infinite systems, see Fig. [2;. The conductance was calculated by performing 
the summation for 100 atoms from each of the first two full layers. As a result we obtained a conductance of 2.58 
Go- Papanikolaou et alM> got a conductance of 3 Go for an infinite Cu wire to be associated with three conducting 
channels within the Landauer approach. For an infinite wire the transmission probability is unity for all states, 
therefore, the conductance is just the number of bands crossing the Fermi level. For the present case of a finite chain, 
the transmission probability is less than unity for all the conducting states. This qualitatively explains the reduced 
conductance with respect to an infinite wire. 

Finally, we studied the dependence of the conductance on the thickness of the nanocontacts. All the investigated 
structures have C± v symmetry and the central layer of the systems is a plane of reflection symmetry. The set-up of the 
structures is summarized in Table [I] Contact refers to a broken contact, while the others have different thicknesses 
from 1 up to 9 Au atoms in the central layer. In Fig. HJthe calculated conductances are displayed as performed by 
taking nearly 100 atoms from each of the first two full layers: I = C — 4, J = G + 4 for the broken contact and 
I = C — 3,J = G + 3 for all the other cases, see Table Q] It can be seen that the conductance is nearly proportional 
to the number of Au atoms in the central layer. This finding can qualitatively be compared with the result of model 
calculations for the conductance of a three-dimensional electron gas through a connective neck as a function of its 
area in the limit of $0 = 90° for the opening angle 31 . In the case of the broken contact, the non-zero conductance 
can again be attributed to tunneling of electrons. 



B. Gold contact with an impurity 

In recent break junction experiments^ remarkable changes of the conductance histograms of nanocontacts formed by 
AuPd alloys have been observed when varying the Pd concentration. Studying the effect of impurities placed into the 
nanocontact are, in that context, at least relevant for dilute alloys. The interesting question is whether the presence of 
impurities can be observed in the measured conductance. For that reason we investigated transition metal impurities 
such as Pd, Fe, and Co placed at various positions of the point contact as shown in Fig. For the notation of the 
impurity positions see Fig. [7| The calculated spin and orbital moments of the magnetic impurities are listed in Table 
ITT1 As usual for magnetic impurities with reduced coordination number—, both for Fe and Co we obtained remarkably 
high spin moments, and in all positions of a Co impurity large orbital moments. In particular, the magnitude of the 
orbital moments is very sensitive to the position of the impurity. This is most obvious in case of Fe, where at positions 
B and C the orbital moment is relatively small, but at position A a surprisingly high value of 0.47 \xb was obtained. 
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The summation over 116 atoms from each of the first two full layers (/ = C — 3, J = C + 3) in Eq. (|16fl has 
been carried out in order to evaluate the conductance. The calculated values are summarized in Table IIIII A Pd 
impurity (independent of position) reduces only little the conductance as compared to a pure Au point contact. This 
qualitatively can be understood from the local density of states (LDOS) of the Pd impurity (calculated by using an 
imaginary part of the energy of (5 = 1 mRyd). In Fig. [SJwe plotted the corresponding LDOS at positions A and C. 
Clearly, the change of the coordination number (8 at position A and 12 at position C), i.e., different hybridization 
between the Pd and Au d bands, results into different widths for the Pd d-like LDOS. In both cases, however, the Pd 
d states are completely filled and no remarkable change in the LDOS at Fermi level (conducting states) happens. 

The case of the magnetic impurities seems to be more interesting. As can be inferred from Table ITTT1 impurities at 
position B change only a little the conductance. Being placed at position A, however, Fe and Co atoms increase the 
conductance by 11 % and 24 %, while at position C they decrease the conductance by 19 % and 27 %, respectively. 
In Ref. ^3 it was found that single Fe, Co (and also Ni) defects in a 2x2 infinite Cu wire decreased the conductance. 
By analyzing the DOS it was concluded that the observed reduction of the conductance is due to a depletion of the s 
states in the minority band. The above situation is very similar to the case of an Fe or Co impurity in position C of 
the point contact considered, even the calculated drop of the conductance (~ -20 % for Fe and ~ -28 % for Co) agrees 
quantitatively well with our present result. Our result, namely, that Fe and Co impurities at position A increase the 
conductance can, however, not be related to the results of Ref. In order to understand this feature we have to 
carefully investigate the LDOS calculated for the point contact. 

In Fig. |5]we plotted the minority <i-like LDOS of the Fe and Co impurities in positions A and C as resolved according 
to the canonical orbitals d x 2_ y 2, d xy , d xz ^ d xy and d 3z 2_ r 2. We have to stress that this kind of partial decomposition, 
usually referred to as the (i,m,s) representation of the LDOS, is not unique within a relativistic formalism, since 
due to the spin-orbit interaction the different spin- and orbital components are mixed. However, due to the large 
spin-splitting of Fe and Co the mixing of the majority and minority spin-states can be neglected. As can be seen from 
Fig. [51 the LDOS of an impurity in position A is much narrower than in position C. This is an obvious consequence 
of the difference in the coordination numbers (8 for position A and 12 for positions C). Thus an impurity in position 
A hybridizes less with the neighboring Au atoms and, as implied by the LDOS, the corresponding d states are fairly 
localized. Also to be seen is a spin-orbit induced splitting of about 8 mRyd (~ 0.1 eV) in the very narrow d x 2_ y 2-d xy 
states of the impurities in position A. The difference of the band filling for the two kind of impurities shows up in a 
clear downward shift of the LDOS of Co with respect to that of Fe. 

In explaining the change of the conductance through the point contact caused by the impurities in positions A and 
C, the s-like DOS at the center site, i.e., at the narrowest section of the contact, is plotted in the top half of Fig. ITU1 
As a comparison the corresponding very flat s-like DOS is shown for a pure Au contact. For contacts with impurities 
this s-like DOS shows a very interesting shape which can indeed be correlated with the corresponding d 3z 2„ r 2-like 
DOS at the impurity site, see bottom half of Fig. ^| Clearly, the center positions and the widths of the d 3z 2_ r 2-like 
DOS peaks and those of the respective (anti-)resonant s-like DOS shapes coincide well with each other. This kind 
of behavior in the DOS resembles to the case studied by Fano for a continuum band and a discrete energy level in 
the presence of configuration interaction (hybridization)^. Apparently, by keeping this analogy, in the point contact 
the s-like states play the role of a continuum and the <i 3z 2_ r 2-like state of the impurity acts as the discrete energy 
level. Since the two kinds of states share the same cylindrical symmetry, interactions between them can occur due to 
backscattering effects. It should be noted that similar resonant line-shapes in the STM I-V characteristics have been 
observed for Kondo impurities at surfaces^*^ and explained theoretically 35 . 

Inspecting Fig. ^3 the enhanced s-like DOS at the Fermi level at the center of the point contact provides a nice 
interpretation to the enhancement of the conductance when the Fe and Co impurity is placed at position A. As the 
peak position of the d 3z 2_ r 2-like states of Fe is shifted upwards by more than 0.01 Ryd with respect to that of Co, the 
corresponding resonance of the s-like states is also shifted and the s-like DOS at the Fermi level is decreased. This is 
also in agreement with the calculated conductances. In the case of impurities at position C, i.e., in a position by a = 
7.63 a.u. away from the center of the contact, the resonant line-shape of the s-like states is reversed in sign, therefore, 
we observe a decreased s-like DOS at the Fermi level, explaining in this case the decreased conductance, see Table 
IIIII Since, however, the s-like DOS for the case of a Co impurity is larger than for an Fe impurity, this simple picture 
cannot account correctly for the opposite relationship we obtained for the corresponding conductances. 
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Summary 



By using a Green's function technique based on the embedding scheme of the multiple scattering theory and the 
Kubo-Greenwood linear response theory as formulated by Baranger and Stoned we studied the conductance of gold 
nanocontacts depending on the contact geometry and transition metal impurities placed at various positions. We 
performed several numerical tests that proved the efficiency of our method. In good agreement with experiments 
and other calculations we obtained a conductance of 1.1 Go for a finite linear chain connecting two semi- infinite Au 
leads. The calculated conductance for a thicker 2x2 wire, 2.58 Go, can be related to a recent result for an infinite 
2x2 chain (3 Gq)^. Also in agreement with quantum mechanical model calculations^! we found a nearly linear 
dependence of the conductance on the "thickness" of the contact. By embedding magnetic transition metal impurities 
into a point contact we found both enhancement and reduction of the conductance depending on the position of the 
impurities. On analyzing the local density of states we concluded that the effect of the impurity is mainly controlled 
by the interaction of the minority <i-like and s-like states giving rise to a resonant line-shape (Fano-resonance) in the 
s-like DOS at the center of contact. We suggest that this line-shape should also be observed in / — V conductance 
characteristics providing thus an "experimental" tool to detect magnetic impurities (even their position) in a noble 
metal point contact. 
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TABLE I: Set-up of various nanocontacts. The table shows the number of Au atoms in the layers as labelled by C, C ± 1, etc. 
in Fig. Contact 1 refers to Fig. 



layer 


Contact 


position 





1 


4 


5 


9 


C±4 


Full 


Full 


Full 


Full 


Full 


C±3 


9 


Full 


Full 


Full 


Full 


C±2 


4 


9 


16 


21 


25 


C±l 


1 


4 


9 


12 


16 


C 





1 


4 


5 


9 



TABLE II: Calculated spin and orbital moments of magnetic impurities placed at different positions in a Au point contact, see 
Fig. 



position 


S z [hb] 






Fe 


Co 


Fe 


Co 


A 


3.36 


2.01 


0.47 


0.38 


B 


3.46 


2.17 


0.04 


0.61 


C 


3.42 


2.14 


0.07 


0.22 



TABLE III: Calculated conductances of a Au point contact with impurities on different positions, see Fig. |7| 



impurity 
position 


Conductance [Go] 


Pd 


Fe 


Co 


A 


2.22 


2.67 


2.97 


B 


2.24 


2.40 


2.26 


C 


2.36 


1.95 


1.75 


pure Au 


2.40 




FIG. 2: Perspective view of some contacts between two fcc(OOl) semi-infinite leads. Only the partially filled layers are shown, 
a) point contact b) slanted linear finite chain c) 2x2 finite chain. 
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FIG. 3: Calculated conductance as a function of the Fermi level broadening, S, for the Au contact shown in Fig. The 
dashed straight line is a linear fit to the values for 6 = 1.5, 2.0, 2.5, and 3.0 mRyd. 




FIG. 4: Conductance versus the number of sites included in the sum in Eq. I Kill for the contact in Fig. The different 
curves show conductances as calculated between different pairs of layers. For a definition of the layer numbering see Fig. 



11 




Number of atoms in the corresponding layers 



FIG. 5: Conductance versus the number of sites included in the sum in Eq. <!16t for the slanted wire shown in Fig. Full 
circles are the results of summing in layers I — C — 4 and J — C + 4 (first full layers) , while squares refer to a summation in 
layers I = C — 2 and J — C + 2 (layers containing 4 Au atoms). 
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FIG. 6: Conductance versus the number of Au atoms in the central layer for the Au contacts described in Tabic U 
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FIG. 7: Impurity positions (light grey spheres) in a Au point contact, see Fig. 




FIG. 8: Local density of states of a Pd impurity in position A (solid line) and in position C (dashed line) of a Au point contact, 
see Fig. Q 
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FIG. 9: Minority-spin orbital-resolved d-like local density of states of Fe and Co impurities in position A (upper panels) and 
in position C (lower panels) of a Au point contact, see Fig. |7] 
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FIG. 10: Top left: minority-spin s-like local density of states at the center site of a Au point contact with an impurity at 
position A, see Fig. [T] (solid line: Co, dashed line: Fe). Top right: the same as before, but with an impuritiy at position C. As a 
comparison, in both figures the corresponding LDOS for the pure Au contact is plotted by dotted lines. The solid vertical lines 
highlight the position of the Fermi energy. Bottom: minority-spin d 3z 2_ r 2 local density of states of the impurities (solid line: 
Co, dashed line: Fe) at positions A (left) and C (right). Vertical dashed lines mark the center positions of the d 3z 2_ r 2-LDOS 
peaks. 



